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Measurements of intensity and polarization of diffuse Galactic synchrotron 
emission as well as starlight polarization reveal power law spectra of fluctuations. 
I"-- We show that these fluctuations can arise from magnetohydrodynamic (MHD) 

turbulence in the Galactic disk and halo. To do so we take into account the 
converging geometry of lines of sight for the observations when the observer is 
^ ■ within the turbulent volume. Assuming that the intensity of turbulence changes 

along the line of sight, we get a reasonable fit to the observed synchrotron data. 
■ As for the spectra of polarized starlight we get a good fit to the observations 

taking into account the fact that the observational sample is biased toward nearby 
O ■ stars. 

Subject headings: cosmic microwave foregrounds-galaxy:structure-MHD- turbulence- 
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1. Introduction 

Attempts to determine the statistics of intensity and polarization fluctuations of cosmic 
microwave background (CMB) renewed interest to the fluctuations of Galactic foreground 
radiation (see Tegmark et al. 2000). Spectra of intensity of synchrotron emission and syn- 
chrotron polarization (see papers in de Oliveira-Costa & Tegmark 1999) as well as starlight 
polarization (Fosalba et al. 2002; henceforth FLPT) have been measured. Those measure- 
ments revealed a range of power-laws, the origin of which has not been addressed as far as 
we know. In Tegmark et al. (2000) there is an allusion that the spectra may be relevant to 
Kolmogorov turbulence, but the issue of how those different spectra may arise has not been 
addressed. 



Interstellar medium is turbulent and Kolmogorov-type spectra were reported on the 
scales from several AU to several kpc (see Armstrong et al. 1995; Lazarian & Pogosyan 
2000; Stanimirovic & Lazarian 2001). Therefore it is natural to think of the turbulence 
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as the origin of the fluctuations of the diffuse foreground radiation. Interstellar medium 
is magnetized with magnetic field making turbulence anisotropic. It may be argued that 
although the spectrum of MHD turbulence exhibits scale-dependent anisotropy if studied 
in the system of reference defined by the local magnetic field (Goldreich & Sridhar 1995; 
Lithwick & Goldreich 2001; Cho & Lazarian 2002), in the observer's system of reference the 
spectrum will show only moderate scale-independent anisotropy. Thus from the observer's 
point of view Kolmogorov description of interstellar turbulence is acceptable in spite of 
the fact that magnetic field is dynamically important and even dominant (see discussion in 
Lazarian & Pogosyan 2000; Cho, Lazarian & Vishniac 2002). 

It is customary for CMB studies to expand the foreground intensity over spherical 
harmonics Y [m , 1(0, (ft) = ^ z m ai rn Yi rn (9 1 0), and write the spectrum in terms of C\ = 
J2m=-i \ a im\ 2 /(2l + 1)- The measurements indicate that angular power spectra (Cj) of the 
Galactic emission follows power law (C t oc l~ a ) (see FLPT and references in §3). 

This paper tests whether the measured spectra are compatible with the theoretical pre- 
dictions of spatial statistics that is expected in the presence of MHD turbulence. Analytical 
studies in this direction were done in Lazarian (1992, 1994, 1995) and here we supplement 
them with numerical simulations of the synthetic spectra. 



In this section we show that, when the angle between the lines of sight is small (i.e. 6 < 
L/d max ), the angular spectrum Ci has the same slope as the 3-dimensional energy spectrum 
of turbulence. Here L is the typical size of the largest energy containing eddies, which is 
sometime called as outer scale of turbulence or energy injection scale, and d max is the distance 
to the farthest eddies (see Fig. la). 

To illustrate the problem consider the observations with lines of sight being parallel. 
The observed intensity is the intensity summed along the line of sight, r z . 



2. Spectra of Fluctuations 



2.1. C\ for small angle limit 




(1) 




Rearranging the order of summation and using f dr z e 



S(k z ), we get 




(3) 
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which means Fourier component of l2D(r x ,r y ) is izD(k Xl k y , 0). 

When the angular size of the observed region (A9 x A9 in radian) on the sky is small, 
Ci is approximately the 'energy' spectrum of hoi^x^y) (Bond & Efstathiou 1987; Hobson 
& Majueijo 1996; Seljak 1997), which means Ci ~ \i3D(k x , k y , 0)| 2 with I ~ k(n/Ad) and 
k — (k 2 + k y ) l l 2 . The analysis of the geometry of crossing lines of sight is more involved, but 
for power-law statistics it follows from Lazarian & Shutenkov (1990) that if \izd\ 2 oc k~ m , 
then the 'energy' spectrum of l2D(r x ,r y ) is also k~ m . Therefore, we have 

Q oc \i3D(k x , ky, 0)| 2 oc r m . (4) 

in the small 9 limit 1 . 

For Kolmogorov turbulence ( |i3z> | 2 oc A; -11 / 3 ), we expect 

Q oc T 11/3 , if 9 < L/d rnax . (5) 

Note that I ~ ir/9. 



2.2. Ci for large angle limit 

Following Lazarian & Shutenkov (1990), we can show that the correlation function for 

9 > Lj d max , 

K{9) = f J dr 1 dr 2 K,{\v l -V2\), 

f°° r^ 2 const 

/ dr r/C(r) / dip ~ — - — , (6) 
Jo Jo & 



sin^ 



where we change variables: (ri,r2) — > (r, ip) is clear from Fig. lb; we accounted for the 
Jacobian of which is rj sin#. We can understand 1/9 behavior qualitatively from Fig. lb. 
When the angle is large, points along of the lines-of-sight near the observer are still correlated. 
These points extend from the observer over the distance oc 1/ sin (9/2) . 



1 In some cases, we infer C; from the observation of the correlation function IC(r) in real space (or 
K{9) on the sky). When the three-dimensional spectrum of turbulence follows a power law (~ k~ m ), 
JC(r) oc /Co — r m ~ 2 , where /Co ~ L m ~ 2 is a constant. However, when the slope of the turbulence spectrum is 
steeper than fc -4 , this relation breaks down and /C(r) oc /Co — r 2 regardless of the turbulence slope. When 
we infer C; from this fC(r), we obtain C; oc Z -4 regardless of the true slope, when the slope is steeper than 
-4. 
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Fig. 1. — (a) Small 9 limit (9 < L/d max ). The fluctuations along the entire length of the lines 
of sight are correlated, (b) Large 6 limit (9 > L/d max ). Only points close to the observer are 
correlated, (c) Angular correlation function for disk (solid line) and for halo (dotted line). 
When 9 is small, K{0) oc constant — 9 5 ^ 3 . When 9 is large, K{9) oc 9~ l . (d) Spectra for disk 
(d max = 10 kpc; turbulence size L = 100 pc) and for halo (d max = 1 kpc; turbulence size 
L = 100 pc). (e) Angular spectrum for Galactic disk. We considered geometry of Galactic 
disk, (f) Angular spectra of degree of polarization. Fixed Distance means all star are at the 
same distance. Random Distance means random sampling of stars according to equation 
(13). Zero point of the y-axis is arbitrary for all graphs. 

In the limit of 9 1 we get the angular power spectrum Ci using Fourier transform: 

Ci ~ J J K(9)e- iie d9 x d9 y 

d9 9J {W)K{9) oc T 1 , (7) 

where 9 = (9 2 x + Oy) 1 / 2 , J is the Bessel function, and we use K(9) oc 
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In summary, for Kolmogorov turbulence, we expect from equations (5) and (7) that 



which means that the power index a of C\ is 2 —1 < a < —11/3. The critical number l cr (~ 
3d max /L) depends on the size of the size of the large turbulent eddies and on the direction of 
the observation. If in the naive model we assume that turbulence is homogeneous along the 
lines of sight and has L ~ 100 pc corresponding to a typical size of the supernova remnant 
for disk with d 

max ~ 10 kpc, we get l cr ~ 300. For the synchrotron halo, d max ~ 1 kpc (see 
Smoot 1999) and we get l cr ~ 30. 



Recent statistical studies of total synchrotron intensity include Haslam all-sky map at 
408 MHz that shows that the power index a is in the range between 2.5 and 3 (Tegmark & 
Efstathiou 1996; Bouchet, Gispert, & Puget 1996). Parkes southern Galactic plane survey at 
2.4 GHz suggests shallower slope: Giardino et al. (2002) obtained a ~ 2.05 after point source 
removal and Baccigalupi et al. (2001) obtained a ~ —0.8 to —2. Using Rhodes/ Hart RAO 
data at 2326 MHz, Giardino et al. (2001) obtained a ~ 2.43 for all-sky data and a ~ 2.92 
for high Galactic latitude regions with |6| > 20°. The rough tendency that follow from these 
data is that a which is close to 2 for the Galactic plane gets steeper (to ~ —3) for higher 
latitudes. In other words — 2 < a < — 3 which differs from naive expectations given by 
equation (8). 

Can the difference be due to non-linear law of synchrotron emissivity? For synchrotron 
radiation, emissivity at a point r, i(r) oc n(e)|-Bj_| 7 , where n(e) is the electron number den- 
sity, B± is the component of magnetic field perpendicular to the line of sight. The index 7 
is approximately 2 for radio synchrotron frequencies (see Smoot 1999). If electrons are uni- 
formly distributed over the scales of magnetic field inhomogeneities, the spectrum of intensity 
reflects the statistics of magnetic field. For small amplitude perturbations (5b/ B <^ 1; this 
is true for scales several times smaller than L when we interpret B as local mean magnetic 
field strength), if 5b has a power- law behavior, the statistics of intensity will have the same 
power-law behavior (see Deshpande et al. 2000). Therefore emissivity non-linearity does not 
account for the difference between the observed spectra and the theoretical expectations. 



2 Note that point sources would result in a ~ 0. 




(8) 



3. Synchrotron Radiation 



- 6- 



To address the problem, we perform simple numerical calculations for galactic disk and 
halo. We obtain C t using the relation 

K(cos9) = -^V(2/ + l)C,P,(cos0), (9) 

I 

c i = \J Pi(cos9)K(cos9) d(cos9). (10) 

We use Gauss-Legendre quadrature integration (see Szapudi et al. 2001 for its application 
to CMB) to obtain C x . We numerically calculate the angular correlation function from 

fT(cos#) = JdnJ dr 2 /C(|ri - r 2 |), (11) 

where |ri — r 2 | = r 2 + r 2 — 2r^ 2 cos9 and assume that the spatial correlation function IC(r) 
follow Kolmogorov statistics: 

w x / /C -r 2 / 3 ifr<L 

M r ) S n t r ( 12 ) 

v ' \ if r > L, v 7 

where /Co ~ £ 2//3 is a constant. Fig. lc and Fig. Id illustrate the agreement of our calculations 
with the theoretical expectations within the naive model of the disk and the halo from the 
previous section. 

To make the spectrum closer to observations we need to consider more realistic models. 
First, synchrotron emission is stronger in spiral arms, and therefore we have more syn- 
chrotron emission coming from the nearby regions. Second, if synchrotron disk component is 
sufficiently thin, then lines of sight are not equivalent and effectively nearby disk component 
contributes more. Indeed, when we observe regions with low Galactic latitude, the effective 
line-of-sight varies with Galactic latitude. 

Suppose we observe a region with b = 10°. Then emission from d — 10 kpc is sub- 
stantially weaker than that from d = 100 pc, because the region with d — 10 kpc is 
10 kpc sin 10° ~ 1.7 kpc above the Galactic plane and, therefore, has weak emission. To 
incorporate this effect, we use the weighting function w(r) = 100/maa;(100, r sin 10°), which 
gives more weight on closer distance. The resulting angular power spectrum (Fig. le) shows 
a slope similar to -2. 

For halo, the simple model predicts that Q oc Z~ n / 3 , but observations provide somewhat 
less steep spectra. Is this discrepancy very significant? The spectrum of magnetic field is 
expected to be shallower than k~ u ^ 3 in the vicinity of the energy injection scale and at the 
vicinity of the magnetic equipartition scale. The observed spectrum also gets shallower if 
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dmax gets larger. For instance Beuermann, Kanbach, & Berkhuijsen (1985) reported the 
existence of thick radio halo that extends to more than ~ 3 kpc near the Sun. Finally, 
filamentary structures and point sources can make the spectrum shallower as well. Further 
research should establish the true reason for the discrepancy. 

4. Galactic Starlight Polarization 

Polarized radiation from dust is an important component of Galactic foreground that 
strongly interferes with intended CMB polarization measurements (see Lazarian & Prunet 
2001). FLPT attempted to predict the spectrum of the polarized foreground from dust and 
obtained C\ ~ Z~ L5 for starlight polarization. This spectrum is different from those discussed 
in the previous sections. To relate this spectrum to the underlying turbulence we should 
account for the following facts: a) the observations are done for the disk component of the 
Galaxy b) the sampled stars are at different distances from the observer with most of the 
stars close-by. 

To deal with this problem we use numerical simulations again. We first generate a 
three (i.e. x,y, and z) components of magnetic field on a two-dimensional plane (4096 x 4096 
grid points representing 10 kpc x 10 kpc), using the following Kolmogorov three-dimensional 
spectrum: E 3D (k) oc k~ u ^ 3 if k > k , where k ~ 1/100 pc. Our results show that how we 
continuously extend the spectrum for k < k does not change our results. 

To simulate the actual distribution of stars within the sample used in FLPT, we scatter 
our emission sources using the following probability distribution function: 

P(r) oc e - r/1 - 5kpc . (13) 

The starlight polarization is due to the difference in absorption cross section of non-spherical 
grains aligned with their longer axes perpendicular to magnetic field (see review by Lazarian 
2000). In numerical calculations we approximate the actual turbulent magnetic field by 
a superposition of the slabs with locally uniform magnetic field in each slap and assume 
that the difference in grain absorption parallel and perpendicular to magnetic field results 
in the 10% difference in the optical depths ry and r_|_ for a slab. We calculate evolution 
of Stocks parameters of the starlight within the slab and use the standard transforms of 
Stocks parameters from one slab to another (Dolginov, Gnedin, & Silantev 1996; see similar 
expressions in Martin 1972). 

We show the result in Fig. If. For comparison, we also calculate the degree of polariza- 
tion assuming all stars are at the same distance of ~ 10 kpc. The result shows that, for a 
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mixture of nearby and faraway stars, the slope steepens and gets very close to the observed 
one, i.e. —1.5. 

Our conclusions may be tested if stars at particular distance only are correlated. If 
those stars are nearby, we would expect to have steeper index and it will become shallow if 
only distant stars are chosen. Choosing only nearby stars we expect to get a steeper index. 
A systematic study of this change can provide an estimate of the energy injection scale L. 
From the point of view of foreground studies, we conclude that a naive use of the polarization 
template produced with the random sample of stars may be misleading. 

5. Summary 

In this paper we have addressed the origin of spatial fluctuations of emissivity and polar- 
ization of Galactic diffuse emission. We have shown that MHD turbulence can qualitatively 
explain observed properties of total synchrotron emission and starlight polarization. The 
variety of measured spatial spectra of synchrotron emission can be explained by the inhomo- 
geneous distribution of emissivity, structure of Galactic disk and halo, and/or various energy 
injection scales. Similarly, MHD turbulence plus inhomogeneous distribution of stars can 
explain the reported scaling of starlight polarization statistics. 

Evidently more systematic studies are required. Those studies will not only give insight 
into how to separate CMB from foregrounds, but also would provide valuable information 
on the structure of interstellar medium and the sources/energy injection scales of interstellar 
turbulence. 

We acknowledge the support of NSF Grant AST-0 125544. 
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